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Abstract 



The case of a hard type II superconductor in the form of strip with elliptic 
cross-section when placed in transverse magnetic field is studied. We approach 
the problem in two steps, both based on the critical-state model. First we 
calculate numerically the penetrated current profiles that ensure complete 
shielding in the interior, without assuming an a priori form for the profiles. 
In the second step we introduce an analytical approximation that asumes that 
the current profiles are ellipses. Expressions linking the sample magnetization 
to the applied field are derived covering the whole range of applied fields. The 
theoretical predictions are tested by the comparison with experimental data 
for the imaginary part of AC susceptibility. 



I. INTRODUCTION 



The successful use of high-T c superconductors in the fabrication of conductors for high 
current applications requires a deep understanding of their response to the magnetic field 
and, in particular, their ac losses behavior. An important case to study is the magnetic 
response of a tape made from hard type II superconductor, with a cross-section of elliptic 
shape, to a magnetic field applied perpendicular to its major axis. This is the configuration 
very often met in practice and includes the common case of a cylindrical wire in perpendic- 
ular field. There were several attempts to study this system, all based in the critical-state 
model Q. Within this model framework, the application of a magnetic field in a supercon- 
ductor results in the penetration of current with a constant density J c . The current profile 
corresponding to a given applied field is distributed in such a way that it keeps the field 
in the internal region unchanged (so, in the initial magnetization curve, the interior field is 
kept zero). Once a method is found to obtain the actual shape for the current profiles -also 
called flux fronts-, all the magnetic properties such as magnetization, ac susceptibility and 
ac losses can be deduced from it. 

A first approach to the problem of a superconducting tape with elliptical cross-section 
assumed that the flux fronts are ellipses with one constant axis (that in the direction of the 
applied field) and the other axis varying in order to match the field change However, 
the assumption of the constant axis was shown to be incorrect in computations of cylindrical 
wires in transverse field f|, spheres and spheroids |5|,|] and also thick strips @], where it 
was shown that field profiles detach from the surface. 

Another important point is the actual shape of the profiles. Ashkin 0] developed a 
numerical technique to calculate current penetration profiles by forcing the field on the 
current penetration profile to vanish. Although the magnetic field in the current-free region 
was zero, as required by critical-state model, with a precision of only around 10 %, the current 
profiles results showed a tendency towards a spindle shape (similar to the profiles shown in 
solid line in Fig. 1). This fact was further pursued by Kuzovlev []8| who demonstrated that, 
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at least for the case of a sphere, the current profiles have indeed a spindle shape, and not a 
more smooth shape such as an ellipsoid. However, Bhagwat and Chaddah solved the case 
of a very thin elliptical tape assuming elliptical flux fronts Q. 

Therefore, the situation on such an important system remains unsolved, and questions 
arise as what the actual shape of current profiles is (elliptical or spindle) and how correct 
and accurate the various approximations presented up to now are. In this work, we will use 
a numerical procedure based on the critical-state model to accurately determine the flux 
fronts that shield the central region of the sample, for any applied field and sample aspect 
ratio. A key feature of our approach is that, different from the above mentioned models, 
we will not assume any a priori shape for the flux fronts. After briefly describing the 
main characteristics of the calculated profiles, we will introduce an analytical model which 
reproduces the features of the actual profiles with enough precision for most needs. We will 
finally compare our results with experimental data measured on an actual superconducting 
tape. 



II. NUMERICAL RESULTS 

Our numerical model is based on minimizing the magnetic energy of the current distri- 
bution after each applied field variation. This approach has been successfully applied to 
describe the experimental features observed in the initial magnetization [T]J as well as the 



whole magnetization loop and levitation force [12] of superconducting cylinders. The details 



of the model can be found in 11.12 



Calculated profiles are shown in Fig. 1 (solid lines) for the cases of superconducting 
tapes of elliptical cross-section with 6/a=0.1 (thin ellipse), 1 (circular wire), and 10 (long 
ellipse), where a and b are the ellipse semiaxes. The external field is applied in the direction 
of the (vertical) b axis. In order to confirm the validity of our approach, we checked that 
the field in the non-penetrated region is zero with a precision of around 0.1 %. As a further 
check of the model, we have found that our calculations for thin elliptical tapes with the 
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aspect ratio smaller or equal than b/a = 0.01 coincide within numerical accuracy with the 
analytical formulas for thin ellipses of Bhagwat and Chaddah ||. This agreement is not 
observed for thicker samples. Another characteristic of the data is that they follow the 



known dependence proportional to cos 9 predicted for round wires in low applied fields [10 



Finally, it is important to remark that the results from our numerical approach have been 
further confirmed by calculations based on the analogous, although independent, approach 
proposed by Brandt |7j. 

The calculations show several interesting features, which allow us to answer the unsolved 
questions posed above. First, all the profiles are spindled-shaped, although when the ellip- 
tical cross-section of the tape is large in the direction of the applied field (case b/a = 10, 
for example) the shape of the profiles resemble more that of an ellipse. Also, in all cases 
the profiles detach from the surface in the b axis. The effect is more clear for thick samples 
than for thin ones. This means that the classical approach of Wilson [tj], which assumed no 
detachment, could work for thin samples but it is not expected to describe accurately the 
situation of a cylindrical wire in perpendicular applied field. Another interesting feature is 
that the spacing between the successive field fronts is rather constant for the thicker ellipses 
but not for the thinner ones. 



III. ANALYTICAL APPROXIMATION 

After having obtained the adequate numerical description of the actual current profiles for 
the superconducting tape with elliptical cross-section, we would like to know if these results 
could be interpolated by an opportune analytical expression. Here the motivation stems 
from the fact that analytical approaches have important practical advantages. Therefore 
we worked out an analytical model that, in spite of its simplicity, reproduces the magnetic 
results of the numerical model with sufficient accuracy. 

The model is based on the assumption that flux fronts are ellipses with semiaxes a^e 
and boE 1 /™, where a and b are the semiaxes defining the strip's cross-section (see Fig. 2). 
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Each ellipse is characterised by the independent variable e, with the values ranging from 1 
(flux front coinciding with the surface) to (flux front collapsed to the centre of the ellipse). 
To generalise the previous approaches, an additional parameter characterising the shape of 
the flux front, n,was introduced . Indeed, the limit of n — > oo corresponds to the constant 
axis model used for a round wire 0, while n = 2 would reproduce the assumptions in ||. 
We will see that the approximation of flux fronts for tapes with different aspect ratios will 
require to adjust n accordingly. 

In the critical state approach we utilise here, the current distribution in an elliptic strip 
with aspect ratio j3 = bo/ao is expressed in polar coordinates as (Fig. 2) 

+j c for r £ < r < r\ f| -n/2 < 9 < n/2 

for r < r e [j r > r\ (1) 
—j c for r £ < r < r x f| n/2 < 8 <3n/2 



where the outer shape of the strip is defined by r\{9) = b / J [3 2 cos 2 (6) + sin (6) and the 



flux front is given as r £ (6) = £& o /V/3 2 cos 2 (0) +e 2 - 2 / n sin 2 (#) . The current distribution (0) 
generates in the ellipse center the magnetic field in the y-direction 

H y( £ ) = cosd ( r i( ) ~ r e(Q)) d0 

= -H p {g{\) - g{e)) (2) 



where H p = 2J c bo/n is the penetration field for a round wire of radius bo, and g{x) is the 
auxiliary function 



x 1/n arctan A / -J=qz - 1 



9(x) = , J x _^ , (3) 



P 2 

,.2-2/r, 



which depends also on parameters (3 and n ||13|| . An important limit for any /3 and n is 



g(0) = 0. The magnetic moment per unit length of the strip with the current distribution 
(pi) extended to length I » a , b in both +z and — z directions is calculated as 



m(e) a ■ f a ° i. x 
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Then, after dividing by the tape cross-section we obtain for the magnetization 

M(s) = = ^(l- e 2 ^ n ) (5) 



7ra b l 3/3 

In the case when the whole section of the strip is saturated with the critical current density, 
e = and one can find the penetration field H s = —H y (0) = H p g(l) and the saturation 
magnetization M s = |M(0)| = 2H p /3f3. 

Let us now describe the magnetization of the elliptical tape. For a zero-field cooled 
sample, in the virgin stage of magnetization, the distribution of currents, which is determined 
by the parameter e , should be such that it shields a field Hq 111 equal to the applied field 
H a . From this condition and Eq. (||) the implicit relation linking H a with Eq in the initial 
magnetization curve is found as 

H?* = H p \g(l) - gfa)] H° hi < H s 

Si = H? M > H s (6) 

with H a = H? hi and % = 0. The initial magnetization M , as a function of £q, is given by 
Eq. (D as 

2# P A 2+1 



M ( eo ) = --g^(l- e r 5 J (7) 

We now analyze the dynamics of flux penetration at the AC field H a = H m cos(ujt) , where 
the maximum field H m corresponds to a e m defined as in Eq. (|6|) for i = m and H^ 1 = H m , 
so that g(e m ) = g{l) — H m /H p if H m < H s , and e m = otherwise. 

The amplitude susceptibility, defined as the ratio of the magnetization and the applied 
field at H a = H m |TJ is 



M (e m ) = 2 (R) 
Xa H a (e m )~ 3(3 g(l) - g(e m ) {) 

At very small amplitudes, the strip is shielded by currents that flow only in a thin surface 
shell and e m — > 1 . The absolute value of the amplitude susceptibility in this limit, denoted 



Xo HH is then 



. . 2 (2n + l)((3 2 - 1) 

3p (n/3 z — l)g{±) — n + 1 

The shielding at low penetrations is a current distribution for which an analytic solution is 

known [[HJ. ^From this solution it follows that the value of Xo is Xo,anai y tic = 1 + 1/(3. The 

correspondence of the value of xo from (|9]) with that calculated according to the analytic 

solution could be then used as a criterion in finding the optimum value of parameter n in 

our model for each value of the parameter (3 as 

2/3 + 3g(l) -5 



(10) 



opt 3(3 2 g(l) - 4(3 + I 
When the applied field is descending (from H m to —H m ) currents in the opposite sense 

begin to enter from surface to inside the tape. Whereas the already present currents are 

kept frozen the new entering currents should shield a field H[ = (H m — H a )/2 P|,P7|. The 

current profile which shields this field is defined by a ex similarly as in Eq. (||) for % = 1. 

So, the magnetization for a given applied field H a (a given H( ht ), considering both new and 

frozen currents, is 

Mi(ei) = M (e m ) - 2M (e 1 ) (11) 

= ^(d +i/n +i-2^ + n 

where S\ is related with the magnetic field as g{s\) = g(l) — (H m — H a )/2H P if H a > H*, and 
E\ — otherwise, and H* = H m — 2H S is the penetration field of the reverse supercurrents 
fl8|j . Notice that if H m < H s , during the reversal stage, reverse currents will never surpass 
the initial ones, S\ > e m > 0. Analogous magnetization expressions can be easily derived 
for the ascending part of the cycle (from — H m to H m ). 

We can now proceed with the calculation of x'' \ the imaginary part of the complex AC 
susceptibility defined as |1| 

X=^T I"™ M 1 (H a )dH a (12) 

7Tjj J—H m 



This quantity allows to calculate the AC loss |p!5 |. The integral (0) is easily determined in 
a numerical way, by inserting the corresponding expressions for the magnetization Mi(H a ) 
for each stage. 



The usefulness of the analytical model described above is that it provides a very con- 
venient approximation to the numerical calculations. In Fig. 1 we show (dashed lines) 
the current profiles corresponding to the numerical ones. The agreement between both is 
very good except in the region closest to the b axis. However, this region is the one that 
contributes less to the magnetic moment of the sample. Significance of neglecting the differ- 
ence between the numerical profiles and the analytical approximation is at best evaluated 
by comparing the results for the magnetization and AC susceptibility. We found for these 
quantities that the results calculated from the analytical model are hardly distinguishable 
in practice from the accurate numerical data. 



IV. EXPERIMENTAL VERIFICATION 

Experimental test of these models were performed on a monocore tape from high-T c 
superconductor Bi 2 Sr2Ca 2 Cu30io in silver matrix. The data are presented in Fig. 3, where 
we show x" as function of the amplitude of the AC field, H m . To compare better the shapes 
of the experimental curves measured at different temperatures with the model proposed 
here, all the curves were normalized to meet in the maximum point of the x" curve. Our 
theoretical results correspond to both the analytical and numerical model (they cannot be 
distinguished in the scale of the picture). It is important to remark that this is a zero- 
parameter fit, since the value of (3 is obtained from the actual sample dimensions and n 
is chosen from f3 after Eq. (0). We see clear distinction between the data obtained at 
superimposed field with respect to those measured without DC field. We explain this by 
the known fact that applying the DC field much larger than the AC field amplitude limits 
the actual magnetic fields to a narrow interval on the j c {B) dependence, approaching the 
assumption of field-independent j c . Thus, only the curves measured with superimposed 
DC field should be compared with our model derived under the assumption of constant j c . 
Indeed, we see that the data registered with superimposed DC field coincide nicely with the 
curve predicted by our model in the region about the maximum. Slight deviations observed 
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at low fields could be attributed to sample imperfections. On the other side, at large AC 
amplitudes - up to 0.015 T used in our experiments - the assumption about constant j c 
does not hold anymore because a wide range of local magnetic fields could lead to quite 
different actual values of j c . This, in our opinion, explains also the huge deviation of data 
measured in zero DC field from our prediction - one can expect under these circumstances a 
significant deformation of flux front shapes due to the dependence of critical current density 
on the magnetic field. Similar conclusion was drawn also in another paper tackling the same 



problem p0[ . The narrowing of the x"(H a ) curve due to j c (B) was predicted for slabs |L9 



thin films [21] and cylinders ||22||. 



V. CONCLUSIONS 

We have numerically calculated the current penetration profiles of a strip of elliptical 
cross-section in a perpendicular applied field, solving some open questions about the shape 
and properties of the profiles. Moreover, we have presented evidence that, for tapes with any 
aspect ratio, assuming that the flux fronts are ellipses where the axis in the field direction 
shrinks as the power 1/n with respect to the shrinking of the perpendicular axis is a good 
approximation, so they provide good basis for a simple calculation of such properties as 
magnetization loops, susceptibilities and AC losses. 
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FIGURES 

FIG. 1. Current profiles for strips with elliptical cross-section of semiaxes b/a = 0.1, 1, and 10, 
corresponding to applied fields H a = 0,0.2,0.4,0.6,0.8, and 1, in units of the penetration field H s 
(from surface inwards). The strip cross-sections have been scaled as circles. Solid lines correspond 
to the numerical method (Sec. [n|), while dashed lines to the analytical approximation (Sec. |III| ). 

FIG. 2. Sketch of the cross section of the elliptical tape. 

FIG. 3. Imaginary part of the AC susceptibility, x", as function of the amplitude of the applied 
ac field, H a . The solid line corresponds to theoretical values, and symbols to experimental data 
obtained from the tape shown in the inset (/? = 0.098), for different temperatures and dc field 
values. 
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